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The original proposal for this grant was to study transition metal triatomics and 
embedded clusters. The interest in this problem arose from two areas. The first 
was an interest in the hydrogen embrittlement problem and the desire to extend 
previous studies of the interaction of H atoms with Fe clusters, which had been 
carried out with explicit inclusion of only the Fe 4s electrons, by including within 
the cluster several Fe atoms with the 4s and 3d electrons present. In addition to 
the materials science applications, the small clusters are also of direct interest, since 
experimental studies of these molecules are being carried out both in matrices and 
in the gas phase. The work on transition metal triatomics was an extension of work 
which had been carried out on transition metal dimers and on Fe clusters under 
an earlier NASA grant (NCC2-148). During the course of this work it became 
important to respond to agency needs in the area of reacting flows. Since the 
author had had considerable experience in computing potential energy surfaces for 
chemical reactions, the research effort was gradually redirected toward potential 
energy surfaces and for the last year of the grant the primary effort has been in 
that area. Thus, this report will describe work in two different areas. The first half 
of the report w'ill concentrate on the transition metal triatomics, while the second 
half of the report will concentrate on chemical reactions. Since much of this work 
has been published this report will be mainly a summary. 

Under NCC2-148, calculations had been carried out on the bonding in the tran- 
sition metal dimers including studies on Sc 2 , Ti 2 , V 2 , Cr 2 , TiV, CU 2 , Y 2 , Nb 2 , and 
M 02 . A summary of this work has been given in an invited paper[l] which was 

T 

written during this grant. An important feature of the bonding in the transition 
metal dimers is that metal-metal d bonding is important for the elements with half 
filled or lower occupancy d shells, but the bonding mainly involves the s electrons 
for those elements with full or nearly full 3d shells (e.g. Cu and Ag). As we shall 
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see these trends also carry over to the transition metal trimers. 

In order to understand the d bonding in the transition metal trimers studies were 
carried out for Ca3, Sc 3 , Scg , and Tig [2] and for Y 3 [3] . The Sc 3 and Y 3 molecules 
were selected for study because matrix ESR experiments had been carried out by 
Weltner and coworkers which established that the ground state of Sc 3 is of doublet 
spin. The latter result indicates that the 3d electrons must be paired into bonding 
orbitals, since three ground state 4s 2 3d 1 Sc atoms lead to a quartet spin if the 3d 
electrons are not paired. From the calculations, it was found that the bonding in 
Ca 3 arises by 4s— >4p promotion and leads to a well depth of about 0.5 eV for R(Ca- 
Ca) = 7.5 ao- For Sc 3 the 4s bonding is similiar to that in Ca 3 , but the 3d electrons 
are also strongly bonding leading to a 2 Ag ground state with a well depth of about 
1.0 eV and R(Sc-Sc) = 5.75 ao- The good 3d bonding orbitals (bonding between all 
three atoms) are 3dag derived from atomic 3d7r" and 3da' x derived from atomic 3d<7, 
while 3d7r' atomic orbitals lead to 3de' orbitals which are bonding between pairs 
of atoms, and the 3dd' and 3dd" derived levels are non-bonding. (Here the atomic 
symmetry is given with respect to an axis connecting the atom to the center of the 
molecule.) Based on the Sc 3 calculations and calculations on Tig , the bonding in 
Vg and Cr^ was also discussed[2]. 

Similar calculations were carried out for Y 3 [3]. Comparison of the wave functions 
for Y 3 and Sc 3 indicates substantial 4d-5p hybridization in Y 3 , but little 3d-4p 
hybridization in Sc 3 . The increased 4d-5p hybridization leads to stabilization of 
4d7r bonding with respect to 4dcr bonding for equilateral triangle Y 3 , and also leads 
to 4d-5p bonding for linear geometries. These effects lead to a different ordering 
of states for equilateral triangle geometries and a smaller excitation energy to the 
linear configuration for Y 3 as compared to Sc 3 . 

The ground state surface and several excited states of Cu 3 were studied at the 
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single configuration SCF plus singles doubles Cl level[4]. Here the bonding mainly 
involves the 4s orbitals but correlation of the 3d electrons leads to significant bond 
contraction effects and 33 electrons were correlated in these studies. These studies 
neglected relativistic effects and also are not size consistent, which presents signif- 
icant problems for CU3 which has three closed shell 3d 10 configurations. In spite 
of these difficulties, which were recognized at the time, this study provided a more 
accurate ground state surface than had previously been available and provided in- 
formation on the excited states which aided in understanding the R2PI spectrum 
of CU3 which had been studied by Morse and Smalley. 

The work on CU3 generated some controversy, since the computed results were 
not in accord with the interpretations of the experimental results which were ac- 
cepted at the time. The upper state in the Morse and Smalley R2PI spectrum 
had been analyzed in terms of a weakly Jahn-Teller distorted upper state of 2 E" 
symmetry, but the calculations favored a 3s Rydberg upper state of 2 A\ symmetry. 
This assignment was supported by the computed upper state geometries, vibra- 
tional frequencies, and oscillator strengths. Since then Morse has reanalyzed his 
spectrum along with new spectral results from the dispersion fluorescence spectrum 
of Rohlfing and Valentini and concludes that the upper state is indeed 2 A\. 

An additional area of controversy was the computed symmetric stretch frequency 
for the ground state of CU3. Theory[4] obtained 224 cm -1 , while Moskovits and 
coworkers observed a progression in 350 cm -1 in the resonance Raman spectrum, 
which they assigned as the symmetric stretch frequency. Since the computations 
neglected relativistic effects and higher excitation effects, it was expected that this 
frequency might be somewhat too low, but by comparison of the computed and 
experimental vibrational frequency for Cu 2 , it was concluded that the symmetric 
stretch frequency of Cu 3 could not be as large as 350 cm -1 . Since then Truhlar and 
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coworkers have combined the computed potential surface for CU3 and the observed 
ground state vibrational levels of CU3 as obtained by Moskovits and by Rohlfing and 
Valentini to obtain an adjusted surface (lower barrier to pseudorotation) which gave 
Jahn-Teller levels in agreement with the levels observed by Rohlfing and Valentini. 
From this analysis it was concluded that the progression in 350 cm -1 observed by 
Moskovits et al. is a progression in Jahn-Teller levels and the symmetric stretch 
frequency is about 270 cm -1 in much better agreement with the results of calcula- 
tions. New calculations by Langhoff, Bauschlicher, Walch, and Laskowski[ 5 ] which 
incorporate the effects of relativity and higher excitations lead to 268 cm -1 in good 
agreement with the new interpretation of experiment. 

Selected portions of the ground state potential energy surfaces of the CU3, Ag3, 
AgCu2, and AuCu2 trimers were studied at the single-reference singles plus doubles 
configuration interaction and coupled pair functional(CPF) levels correlating 33 
electrons[6]. The CPF method incorporates the effects of higher excitations, while 
relativistic effects are incorporated using relativistic effective core potentials. 

The bonding in CU3, Ag3, AgCu2, and AuCu2 molecules is found to arise basically 
from the (n+l)s 1 nd 10 ground states of the atoms as had been found for Cus[ 4 , 5 ]. 
The ground state surface of Ag3 is found to be similar to that of CU3. For AgCu2 
and AuCu2, since the atoms are inequivalent, there are several possible isomers. 
Here the lowest energy structure has the unique atom in the center. Since the Cu- 
Cu bond length is shorter than the Ag-Cu or Au-Cu bond length, these molecules 
have acute angle structures which leads to 2 Ai ground states. For AuCu2 the acute 
angle C 2v structure is a true minimum, but for AgCu2 the potential is very flat with 
respect to distortions leading to C e symmetry. These differences between AgCu2 
and AuCu 2 are related to the larger bond strength for an AuCu bond as compared 
to an AgCu bond. 
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For the Cu atom the 4s 2 3d 9 state is at 1.49 eV, while for Ag the 5s 2 4d 9 state is 
at 3.74 eV, and for Au the 6s 2 5d 9 state is at 1.74 eV. For CU 3 there is significant 
admixture of 4s 2 3d 9 character leading to a 3d population of 29.47 at the CPF level. 
As has also been found for the transition metal hydrides, SCF does not allow for a 
proper mixing of atomic configurations and this defect is not corrected completely 
by singles and doubles Cl. Higher excitations are therefore important for CU 3 and 
CPF leads to significant differences as compared to SDCI. Ag 3 on the other hand 
is mostly 5s J 4d 10 because of the larger excitation energy to the excited state. This 
leads to a larger d population, 29.76 at the CPF level, and higher excitations are 
less important for Ag 3 than for CU 3 . Thus SDCI and CPF give similar results for 
this molecule. For Au on the other hand, the excitation energy decreases again and 
this is reflected in a low d population in AuCu 2 of 29.39 electrons at the CPF level. 
From this result we predict that higher excitations would be important for AU 3 . 

Calculations were also carried out for the Ni 3 molecule. Here the the results of 
experiment and of existing calculations were still ambiguous. Resonance Raman 
experiments had been interpreted in terms of a bent molecule with an apical angle 
of about 100 degrees. This structure was thought to arise from a A u state of 
unspecified multiplicity via either matrix effects or Renner- Teller distortion. 

If only 4s 1 3d 9 states of Ni atom are involved in Ni 3 , we expect that Ni 3 will be 
a near equilateral triangle Jahn-Teller distorted system much like CU 3 . However, 
there is another possibility which was first suggested by Walch and Goddard. For 
the linear configuration if the center Ni atom is promoted to the nearly degenerate 
4s 2 3d 8 state then it is possible to form two bonds to adjacent 4s J 3d 9 Ni atoms. 
We expect the latter state to be low-lying and the calculations show that for linear 
geometries the 4s 2 3d 8 center atom state is slightly (less than 0.1 eV) below the 
all 4s 1 3d 9 structure for near equilateral triangle geometries. From comparison of 
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the computed and observed vibrational frequencies it is concluded that the bent 
structure is the ground state in the matrix, but the linear structure may be the 
ground state in the gas phase. Note that this conclusion is different from the 
speculation above that the lowest state of Ni3 arises by distortion from a linear 
structure. 

An important problem in this area has been the difficulty in doing spectroscopy 
on these molecules. These difficulties arise because the small binding energies of 
these molecules lead to predisociation in most of the excited states. As for CU3, Ni3 
is found to have a 3s Rydberg upper state. Here the transition is predicted to be in 
the range 2. 7-3. 5 eV which is consistent with the work of Gole who finds a strong 
transition in Ni3 at 2.8 eV. In addition, Gole has measured the symmetric stretch 
and bending frequencies of the ground state of Ni3 and these frequencies are found 
to be in good agreement with the calculated values for the lowest computed state 
of Ni3. As an extension of this work calculations were carried out of the the excited 
states of the Ag3 molecule[8]. 

The work on transition metal trimers[9] and transition metal dimers and other 
systems[lO] was presented at a NATO workshop on transition metal bonding. 

In collaboration with other members of the research group calculations were car- 
ried out for dipole moment of NiH[ll]. This work provided experimental conforma- 
tion of the importance of admixture of atomic states in transition metal bonding. 
These effects had been discussed for example for CU3 above. NiH provides a more 
dramatic example in that the ground state of the molecule has about equal parts 
of 4s 2 3d 8 and 4s 1 3d 9 atomic character. It turns out that the dipole moment is a 
very sensitive function of the mixture of atomic states and the experimentally de- 
termined dipole moment by Field and coworkers confirms the mixed state model of 
the bonding in this molecule. These studies were later extended to other hydrides 
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of the first transition row [12]. 

During the second half of the funding period much of the research effort was 
redirected toward the calculation of potential energy surfaces for chemical reac- 
tions. This new emphasis reflected the importance of being able to compute, from 
first principles, reaction rate constants for applications related to the design of 
the National Aerospace Plane (NASP) and Aeroassisted Orbital Transfer Vehicles 
(AOTV’s). Two groups of reactions were considered. 

The first set of reactions is a subset of the reactions important in H 2 combustion. 
These reactions are important in the design of the Supersonic Combustion RAM 
(SCRAM) jet engine for the NASP. Among reactions which are being studied here 
are: 


H + 0 2 -> H0 2 -» OH + O 

(1) 

0( 3 P) +tf 2 — ► OH + H 

(2) 

OH + H 2 -> H 2 0 + H 

(3) 


The combustion of H 2 in the SCRAM jet engine involves three steps. Initiation 
requires production of H or O atoms. Radical propagation then occurs for H atoms 
via reaction (1) or for O atoms via reaction (2). OH radicals produced by reactions 
(1) or (2) react further via reaction (3). These reactions lead to the production of 
radicals in the combustor region of the SCRAM jet engine. The third step, which 
occurs in the nozzle region of the SCRAM jet engine, is radical recombination. 
Among reactions which are important here are: 


OH + H +M -> H 2 0 + M 

(4) 

H + H + M -> H 2 + M 

(5) 
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Since these reactions evolve a considerable amount of heat, the rates are critical to 
computation of the net thrust production of the engine. The rates of these recombi- 
nation reactions are slow compared to the residence time in the nozzle region, and 
the rates are not well knowm; thus, the computated rates for these reactions are 
crucial to the design of the SCRAM jet engine. 

Calculations have been completed for reaction (2) and are currently underway 
for reaction (l). The work on reaction (2) is described in detail elsewhere[l3]. The 
calculations on reaction (2) focused on the problem of obtaining an accurate barrier 
height. Previous work by Walch et al. on reaction (2) using the POL-CI method 
had given a barrier height of 12.5 kcal/mole, but attempts to improve upon this 
result by simple multi-reference singles and doubles Cl methods lead to a barrier 
height of 16.0 kcal/mole. From comparison of computed and experimental thermal 
rates using the POL-CI surface, it was concluded that the 12.5 kcal/mole POL-CI 
barrier height is accurate. Similar results are seen for F + H 2 where first order Cl 
gave a barrier height of 1.66 kcal/mole, but multi-reference singles and doubles Cl 
leads to barrier heights greater than 3.2 kcal/mole. 

The problem of computing the barrier height for these reactions is found to be 
related to the problem of computing the electron affinity of the O and F atoms. 
Ionic character is important in the saddle point and products regions but not in 
the reactants regions for these reactions. In these calculations, the convergence of 
the barrier height was studied as a function of expanding the active space in the 
CASSCF calculation as well as expanding the basis set. It was found that inclusion 
of 2p — ► 2p' terms, which improve the description of the electron affinity of the O 
and F atoms, lowers the barrier height for the O + H 2 and F + H 2 reactions. It 
was also found that the good barrier height obtained using POL-CI results from a 
cancellation of opposing effects; the neglect of angular correlation in the POL-CI 
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favors the saddle point region over reactants while the neglect of ionic terms in the 
POL-CI favors the reactants over the saddle point region. The calculation on O + 
H 2 serves as an important benchmark for reactions involving O atoms. 

For reaction (l) multi-reference singles and doubles Cl calculations have been 
carried out along the minimum energy path for the addition of H atom to O 2 
determined in the POL-CI studies of Dunning et al. The barrier to H atom addition 
is found to be less than 0.3 kcal/mole, in agreement with the earlier POL-CI work. 
Calculations are also being carried out for the process of dissociating HO 2 to OH 
-1- O. In collaboration with C. F. Melius and C. Rohlfing of Sandia laboratories it 
is planned to generate a global potential surface for this reaction. 

The second set of reactions which are being studied are the reactions: 


0 + N 2 

— » NO + N 

(6) 

n + o 2 

— ► NO + O 

(7) 


These reactions are important in high temperature air chemistry. Two areas where 
these reactions are expected to be important are i) the computation of the compo- 
sition of the gas at the inlet to the SCRAM jet for the NASP and ii) the chemistry 
occuring in the shock layer in front of the AOTV. The details of this work are 
described elsewhere[l4]. Calculations for these reactions are essentially complete. 
There has been some difficulty in obtaining a satisfactory fit to the computed points. 
An approximate fit has been obtained by M. Pattengill and R. L. Jaffe and trajec- 
tory calculations on this surface give good agreement with experiment in the high 
temperature region. This surface will subsequently be used to compute rates under 
non-equilibrium conditions such as occur in AOTV applications. 

In addition to the work described here, several collaborative projects have been 
undertaken with other members of the Ames computational chemistry group. One 
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such project was the calculation of the potential curves for the low-lying 3 II U and 
3 £j states of A^IS]- Another related project involved computing the electronic 
structure of the lowest few states of linear and T shaped CuAl 2 [l 6 ]. 

Finally, the POL-CI surface for reaction(2), which the author generated while 
working for the Department of Energy is still being used in dynamics studies[17]. 
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